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Abstract 
CJ . 

O '< We analyze the large-spin asymptotics of a class of spin-network wavefunctions of Eu- 

^ ' clidean Loop Quantum Gravity, which corresponds to a flat spacetime. A wavefunction from 

OX). this class can be represented as a sum over the spins of an amplitude for a spin network 

whose graph is a composition of the the wavefunction spin network graph with the dual 
one-complex graph and the tetrahedron graphs for a triangulation of the spatial 3-manifold. 
^ ' This spin-network amplitude can be represented as a product of 6j symbols, which is then 

used to find the large-spin asymptotics of the wavefunction. By using the Laplace method 
we show that the large-spin asymptotics is given by a sum of Gaussian functions. However, 
these Gaussian functions are not of the type which gives the correct graviton propagator. 
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1 Introduction 

Loop Quantum Gravity (LQG) is a theory of non-perturbative and background-independent 
quantization of GR, see [1]. It is based on the canonical quantization method, and instead of 
spatial metric, a spatial connection is used as the configuration space variable. Consequently, 
the Hilbert space of physical states is spanned by the spin network states \T), where F is a 
closed SU(2) spin network. The graph of F is a combinatorial graph, i.e. it is a homotopy class, 
because the spin-network states are diffeomorphism invariant. A physical state |^) is given as 
a linear combination of the spin-network states such that it satisfies the quantum Hamiltonian 
constraint Til^) = 0, where Ti is the Hamiltonian constraint operator. 

Solving the Hamiltonian constraint is a difficult problem, and various strategies have been 
developed over the years. In particular, one can consider the quantum Hamiltonian constraint 
in the Ashtekar connection representation, where H becomes a polynomial in the functional 
derivatives. Then any functional '^{A) with a support on flat connections is a solution when 
the cosmological constant is zero, see [2] . The Ashtekar connection is complex in the Minkowski 
signature case so that the resolution of the identity is given by 



V{ReA)V{ImA)\A){A\= I, (1) 

since the operator A is similar to the annihilation operator for the harmonic oscillator. The 
resolution of the identity (1) is an obstacle to construct the loop transform from ^{A) to the spin 
network wavefunction ^(F). However, in the Euclidean signature case, the Ashtekar connection 
is real, so that one obtains the usual expression 

fvA\A){A\ =1. 

Consequently 

{r\^) = JvAWriA)^{A), (2) 

where IVr(^) = (^|r) is the generalization of the trace of the holonomy along a curve to a spin 
network, see [1]. 

The path integral (2) can be rigorously defined by using a three-dimensional spin-foam 

state sum for the quantum group Uq{su{2)) where q = exp [-jS^j is a root of unity [2, 3]. 

When ^{A) = const ■ S{F) then (2) becomes an invariant ^^(F), which is proportional to the 
Witten-Reshetikhin-Turaev invariant for a spin network F embedded in a compact 3-manifold 
S representing a spatial slice of the spacetime [4] . This means that the state 

|*,) = J;M/,(F)|F) (3) 

r 

corresponds to the Kodama wavefunction 

*fc(^) = exp (i— / Tr(AAdA + -AAAAA]] . 

Therefore the effect of using a quantum SU{2) group is that the flat-connection state, represented 
by the wavefunction ^{A), is deformed into the Kodama state (3). This state is a physical state 
for quantum GR with the cosmological constant A/, proportional to 1/k, where k is an integer. 
When "^{A) = exp (i J^ Tr{E"^Am)) ^{F), it can be argued that the corresponding quantum 
group state \^k) is a physical state for quantum GR with the cosmological constant A^ such 
that the triads take the values Em, see [3]. Hence the spin-network wavefunction ^fc(F, E) gives 



a physical state describing a spacetime whose spatial metric is given by Qmn = Tr{EmEn)- In 
particular, one can choose the flat triads and hence \^k) can be considered as a Euclidean analog 
of the vacuum state for a De Sitter spacetime. Since A^ — )■ as /c — )■ oo, then \^k) for large k 
can be considered as a good approximation for the Euclidean flat-spacetime vacuum. 

This result is very useful for the problem of finding a semiclassical limit of LQG, since it can 
be shown that the graviton propagator in LQG will have the correct large-distance asymptotics 
if the vacuum wavefunction 'I'(r, E) has a certain kind of Gaussian asymptotics when the spins 
of r are large [5], see also [6, 7] for a related approach. 

As shown in [5] , if the large-spin asymptotics is given by 

M/(r,io)«iV(r,jo)exp|-J]^M£)(j,_j,)(j,_jo)) , (4) 

where ji is a spin of an edge /, jo is a parameter associated to the flat spatial metric and C{T) 
are jo-independent and positive-definite matrices, then the corresponding graviton propagator 
will have the correct large-distance asymptotics. 

In order to analyze the large-spin asymptotics of ^fc(r,jo), we will represent it as a sum 
over spins of products of quantum dimensions and quantum 6j symbols. When k is very large, 
we will approximate ^^ by the corresponding Ponzano-Regge expression, i.e. we will replace 
the quantum dimensions and the quantum 6j symbols in the expression for ^^ with the cor- 
responding classical evaluations. The rationale for this is that ^I'fc(^)io) was constructed as a 
regularization of an expression for the zero-cosmological constant ^(r,jo), which was given by 
the Ponzano-Regge state-sum without the spin cut-off, see [3]. If a spin cut-off k/2 is introduced, 
one obtains another regularization of ^(r,jo), such that the quantum group is not used. This 
was the regularization originally used by Ponzano and Regge. Furthermore, the sums 2^,0 ^'^ 
^fcir^io) will be replaced by the integrals /„ , since k is very large. Consequently 



^(i,Jo)« / rxf{x,j,Jo). (5) 

Jd 

We will split the integration region D as D = I?_ U Z)+ where I?+ is the region where all 
the spins are large (xj > j'o), so that ^ « \1'_ + ^+. Consequently the 6j symbols from \I'_|_ 
can be approximated by the Ponzano-Regge formula, while the 6j symbols from ^_ can be 
approximated by the asymptotic formulas for 5, 4 and 3 large spins, which were also found by 
Ponzano and Regge, see [9]. We will study in detail the asymptotics of ^+, while the analysis 
of ^_ will be only sketched since it is very similar to the analysis of ^-|- and it can be shown 
that the ^_ asymptotic contribution is subleading to that of ^+. 

We will use the Laplace method to find the asymptotics of (5) for large j and jq. In order 
to do this we will first approximate / as a sum of exponentials, see section 3. Then it is not 
difficult to show that 



n,± 

where -B„ e C„ are constant (jo-independent) matrices. 

The expression (6) will give the desired asymptotics if B^ = for all n and all spin networks 
and C^ 7^ for some n of some spin network. In the following sections we will show that either 
-B+ 7^ or C^ = for all n and all spin networks. We will also show that C~ = for all n and 
all spin networks, so that we will prove that ^{j,Jo) does not have the desired asymptotics. 

The paper is organized in the following way. In section 2 we briefly explain the construction of 
the relevant spin-network wavefunctions from [2, 3, 4]. In section 3 we outline the main procedure 



of analyzing the wavefunction asymptotics, and introduce the notation. Section 4 deals with 
the detailed analysis of the large-spin asymptotics of the ^+ part of the wavefunction. The 
^-1- is written as an integral of an exponential function, which is suitable for the stationary- 
point approximation. The integration is then performed in section 5 and the result is a sum 
of Gaussian functions of the form similar to (4), but with a more general matrix coefficient in 
the exponent, denoted as S. In section 6 we analyze this coefficient, and prove that it never 
has the form (4). In order to demonstrate and verify this result further, in section 7 the matrix 
S is explicitly computed for two simple spin-networks, which are a loop spin network and a 
theta spin network. In section 8 we discuss the large-spin asymptotics of the ^_ part of the 
wavefunction. We present our conclusions in 9, while in the Appendix we give all the necessary 
formulas and prove a matrix theorem which determines the asymptotics of the S matrix. 

2 Physical spin-network wavefunctions 

Let S be a compact 3-manifold, and let T = {7, j/, tv} be a spin network embedded in S, where 
7 is the spin-network graph, ji denote the edge spins and l^ denote the vertex intertwiners. 
Let A(S) be a triangulation of S adapted to 7 in the following way: let H he a handle-body 
obtained by thickening of the dual one-complex A^(S). The graph 7 is embedded in H such that 
each vertex v^ is placed in a different 3-handle of H and each edge l^ runs through appropriate 
one-handles of H, see Fig 1. Let Lh be the Chain-Mail link associated with H and let L/ be 
a set of loops associated to the edges of 7, such that Li is embedded in the one-handle of H 
associated to an edge / and L; is linked with the meridian loop for that one-handle, see Fig 2. 




Figure 1. 




Figure 2. 

Let : S — 7- S*^ be a smooth map from S to a 3-sphere, and let us color the link Lh with 
the ri elements. The ^ element is a linear combination of colors given by 

fc/2 

^ = ^dimqjC(i), 

j=Q 

where C{j) denotes the color (spin j) associated to a loop of Lh- We also introduce 

fc/2 
3=0 

Let us color the L/ loops with Uq{su{2)) irreps A/. We will denote the quantum group 
evaluation of the colored link Lh ULiU---UL„U7 embedded in S^ as 



{Lh ULiU---U Ln, U 7, f^"^+"^ , i, i, A) , 



(7) 



where n2 is the number of dual edges (triangles) and 71,3 is the number of dual vertices (tetrahe- 
drons) . 

From the properties of the Q element, see [4], it follows that the evaluation (7) can be 
expressed as 

j'y f 

where {A^ • (Tet)^^ ■ ^,j\L\j,i,X\ is the spin network which is obtained by removing ^2 ^- 
elements from n2 1-handle meridians. 
Then 

^u{T,E)= Y, llK>^i,E,)lldim,j'f{Al.iTetr-^,j',i',j,L,X), (8) 



where 



/i(A,i?) 



1 

dim A 



SU{2) 



dgf{g,E)x'^^Hg). 



and xa is the trace of the A-representation matrix of g. The function f(g, E) is determined by 
the choice of the flat-connection wavefunctional given by 



^{A) = e^v{i j Tr{E"'A„,)\ 



^{A)5{F) 



This '^{A) solves the Hamiltonian constraint for A = 0, but if V'(^) 7^ 1; then ^{A) does not 
have the sharp values for the tetrads, which is the LQG equivalent of replacing a plane-wave 
with a coherent state. 

One can introduce the background spins j^ associated with the triads E, via the relation 
\Ei\ = jfL'^, where 

Ei= [ i?'"e™„pdx"A(ix^ 

A; is the triangle dual to a dual edge / and L is the Planck length. In the flat-triad case, one 
can assume that \Ei\ = const, and therefore ji = jo, so that ^fc(r, E) = ^^(r, jo). We will then 
choose 

"<^' = ^XTr- <') 

in order to mimic the Rovelli ansatz for the wavefunction [6]. We will then examine the asymp- 
totics of (8) for large spins j and Jq. 

We will express the evaluation of the A* • (Tet)"* • 7 spin network as a sum of products of 
quantum 6j symbols, since this will facilitate our analysis of the large-spin asymptotics. This can 
be done because the evaluation of an arbitrary three-valent spin network T can be represented 
as a sum of products of 6j symbols, see [8, 4]. 

Let r' be a projection of F onto an S'^. The graph 7' will divide the sphere into disjoint discs. 
Color the discs with SU{2) irreps ai, ...,«« and write the corresponding Turaev shadow-world 
evaluation w{r,a), which is given as a product of 6j symbols. Then 

n 

(r) oc N^ 1 [ dimg Ui w(r, a) . (10) 

Alternatively, (10) is the quantum group evaluation of the link formed by 7' and the Chain-Mail 
link for a two-dimensional handle-body which is a thickening of 7' [4]. 

By using (8) and (10) we obtain the following expression for the spin-network wavefunction 

*fc(i,jo)= E (n^-n^"n'"(^)ni6i4) , (n) 

where 

• spins j label the graph 7, 

• spins a label the faces of A*, 

• spins A label the edges of AJ, 

• spins a label the disjoint discs of the projected graph A^ • (Tet)"^ • 7, 

• {6ji;} denotes the 6j-symbol associated to a vertex v of the projected graph A^ • (Tet)'^^ -7, 

• dx is the quantum dimension of the representation x. 

• L denote the intertwiners for the graph A* • (Tet)"'^ ■ 7. 

6 



For example, when T = {7,j} where 7 is a loop embedded in S = S^, we can triangulate 
S"^ with two tetrahedrons such that A^ is a theta-four graph 64. The corresponding 3-valent 
64 • {Tet)'^ • 7 spin network is given in Fig 3. in the Appendix. 

We will also use the Latin indices M, N,p, q, r, s to denote any of the j, a, a, A and l indices, 
and 

j ej = {l,...,J}, a£A = {J + l,...,J + A}, 

a,iey = {J + A + l,...J + A + T}, 

X£C = {J + A + T + 1,...,J + A + T + L}, 

N €j\f = {!,.. .,n}, n = J + A + T + L. 

The corresponding sums and products will be over the whole domain appropriate for each type 
of index, unless otherwise noted. We will often have the set of all values of the indices except 
for the j-indices, so we denote it as 

i = Ar/j. 

It has uj = Q — J elements. 

The 6j symbols {6jv} are functions of all spins from a set J. The index v will also be treated 
as a multi-index for spins, in the sense that the 6j symbol 

Xi X2 X3 
X4 X5 Xq 

is enumerated by the index v = (1, 2, 3, 4, 5, 6), representing the ordered 6-tuple of spins which 
appear in that 6j symbol. This will be useful in limiting the domain of indices to a particular 
6j symbol, noted as M, N £ v. 

The function /u(A) can be chosen arbitrarily, but as we explained in section 2, we will choose 
the Gaussian function (9). In this way a spin scale, jo, is conveniently introduced in the wave- 
function. 

3 Preliminary analysis 

As we explained in the introduction, the quantum group expression (11) will be replaced by the 
corresponding classical group expression, where the spins will have a cut-off given by k/2. Since 
we are interested in the case where k is a very large number, then the finite sums ^q in (8) 
can be approximated by the integrals /q , so that 

^ij,h)^ f d^xfix,j,Jo), (12) 

Jd 

where 

/(x,i,io)=n(2«+i)n(2«+i)n^^xTrni6i4- 

a a X V 

The domain D is a subset of {k/2) and D is determined by the triangle conditions for the spins 
coming from the 6j symbols. 

Let us split the integration region D as D = I?_ U -D+ where D^ is the region where all the 
spins are large (x > jo for every x), and D^ = D\ Z?+. Consequently 

^(i,Jo)~/ d^xf{x,j,jo)+ [ d'^xf{x,j,jo) = ^-U,jo) + ^+{JJo). (13) 

JD- Jd+ 



The 6j symbols from \I/_|_ can be approximated by the PR formula, while the 6j symbols from 
\E'_ can be approximated by the asymptotic formulas for 5, 4 and 3 large spins. Namely, there 
will be a certain number of tetrahedrons which contain the large spins j, so that each of these 
tetrahedrons will have at least two other large spins, due to the triangle inequalities. Each of 
these large spins appear in other tetrahedrons, which will force another spins to be large, and 
so on. In the end there will be a substantial number of large spins different from j, but some of 
the internal spins can still remain small^. Consequently 



J D± 

where f± are the corresponding approximations for / in D± regions. 

We will use the Laplace method to find the asymptotics of (12) for large j and jq. In order 
to do this we will first approximate f± as sums of exponentials. Namely, if Xn are the stationary 
points of /(x) then 

/(^,i,Jo)« j;ene-^"("'^'^°), (14) 

n 

where Sn{x,j,jo) « | In |/(a;, j, jo)|| in the vicinity of x„ and e„ = ±1 depending on whether x„ 
is a minimum or a maximum. Consequently 

*(i,Jo) « Ve„ / d^xe-^"(-'^'^°) = Ten In, 

and one can apply the Laplace method to evaluate the integrals /„. This gives 

^ ~ J^ en J^ K(io) e~^0--in)^s;{io)(i-i*) ^ (15) 

n x*,j* 

where x* , j* are the stationary points of S and 

where S"^-, S"* • and S*^ denote the Hessian matrices in the respective stationary points. 

In order to find the matrix functions 5'*(jo) we will use the scaling properties of Sn when 
the spins x,j and j'o are scaled. Namely, let us assume that Sn satisfy 

Sn{Ax,Aj,Ajo) = A^ [Rn{x,j,jo) + 0(1/ A)] . (16) 

The scaling (16) will be consistent with the approximation (15) if the stationary points have the 
form 

X* = /i^io + /uS + O(io~') , f = ^*ih + ^0 + 0{Jo') . (17) 

The form (17) of the stationary points, together with the scaling (16) and the approximation 
(15) imply 

S* = B* + C*j^'+0{j^^), (18) 

where B* and C* are constant (independent of jo) matrices. If Sn are such that z^* = 1, which 
will be imposed by the choice (9), then we will obtain the asymptotics of the same type as (4). 
However, in order to make a final comparison we need to calculate the matrices B* and C*. 



^In the case of the loop spin network, a numerical investigation has given 11.000 configurations with small 
spins. The maximal number of small spins was 9 out of 22 spins. 



4 ^+ integral 

Let us now make a more detailed analysis of the ^+ integral. In order to calculate the matrices 
A* and B* it will be convenient to introduce the scaling parameter A into the integral ^-|- 
through the following change of variables 

1 . 1 . 1 . 

x+- = Ay, j + - = Ayj, jo + - = Ayo . (19) 

Since x = 0(A) and y = 0(1), the integration domain D^^- is then transformed into D'j^ which 
is of 0(1). A is essentially the same as jo, since yo = 0(1) and one can choose yo = 1. The 
Jacobian of the transformation is A'^, while the factors 2a + 1 become 2 Ay a- The spins j and 
jo are changed into new variables yj i yo, by using the formula (19). The exponents e"*-^"-'"' 
become e~ (y>-~yo) ^ so that we obtain 

^+(j,io)«2^A^+- / d^y WvaWvc^Xl \{{(Sj,{Ay)}. (20) 

•^^V a a \ y^ V 

Since every 6j symbol in (20) can now be approximated by the PR formula (57), we obtain 



where 



s,{A,y) = AY,ysesAy) + ^ + \FDLiy) + o(j^] . (22) 



The 0(A) term represents the Regge action, while the explicit form of the complicated 0(1/A) 
term can be found in [10]. 

In order to apply Laplace's method, it is vital to rewrite the integrand as an exponential 
function, with the multiplicative factor A in the exponent, and to determine the positions of 
the extremal points. The main problem is that the integrand is a product of cosine functions, 
which cannot be easily cast into an exponential form. This problem can be solved by using the 
approximation formula (65) derived in Appendix C: 

cosx« ^-^^Y.^-lfe-'^^''-'"'^'. (23) 



??4(0,e- 2 



pez 



This approximation has two main advantages. First, we avoid having to deal with complex- 
valued exponents, which would have been inevitable if we had employed the formula cosx = 
(e*^ + e~*^)/2 and the corresponding stationary-phase method. Using the stationary-phase 
method would make the asymptotic analysis more complicated, because the corresponding sta- 
tionary points will have the coordinates which are complex numbers. Second, we can "capture" 
the neighborhood of all extremal points at once, including those far away — when p ~ 0(A) — 
in a manifest manner. Specifically, every maximum and minimum at infinity can be labeled as 
p = Am + n, where m,n £ 7j are of the order 0(1). 

By applying (23) to (21), we have 



1 
t?4(0,e'^ 






\Arn,„H-n„ -i(5t,-Am„7r-?i„7r)2 



cos5„(A,y) = —^^ Y, (-l)^^'^''+"''e-2^'^"^''""^-""^^ (24) 



2 



mv,nv(i£ 
mt,,n„<CA 



Here it is crucial to note that the approximation is valid iff the exponent goes to zero, which will 
happen in the vicinity of extremal points t/q. We calculate the exact positions of these points 
by using the ansatz 

y*N = A. + ^ + ^ + 0(l,]. (25) 



A A2 



A3 



Substituting this into the exponent and using (22), we obtain the necessary conditions^ for the 
extremal points y^ , in the form of two systems of equations for coefficients An and Bj^ 



'^As6s,viA) = rUyTT, 



(26) 



(27) 



s(^v 



Of course, these equations are not a sufficient condition to determine the extremal points, simply 
because we are yet to discuss the full integrand in (21). The additional missing equations will 
be determined later. 

Returning now to (24), we will use (22) to expand the exponent in (24) in powers of 1/A in 
order to extract the leading A^ term, as needed for the saddle-point method 



-A^ 



+ 



m 



m^TT + 



\.sGv 



He 



ysOs,v - mvTT 



IT 



n^TT + 



\.S&V 



1 /TT 



2A2 V4 



n„TT 



+ A2 ( IZ y^^^'^ ~ "^"^^ \Fdl + iy-^j 



This can be done for all cosine functions in (21), so in the end we obtain 



*+(i,io) 



2^A^+^-^ 



^/l2^T?4(0,e"- 



V 



E 






y (_l)E„(A™.+n.) / d'^ye^'^(^'J'). (28) 



This expression can be explicitly integrated via the Laplace method, term by term. 
Here the phase has the general form 

5(A,y) ^ So{y) + jSiiy) + ^52(y) + O (^ 

and we have explicitly 

Si{y) = - X^ Xl ysOs,v - m^T^ M ~ ~ "-^'^ 

S2{y) = X^lny^--^ InVy + (- - n^TTJ + 2 ^^ ys9s,v - m^ir Fpi 

a V [_ \s£v / 



(29) 



]ysOs,v -rriyTT 



(30) 



^Here we can also note the following detail. In principle, we could have written the integer p from equation (23) 
in the more general form, which includes some higher power of A, like p = A Z + Am + n. In that case we would 
obtain a consistency condition Z = in addition to equations (26) and (27). Therefore, our choice p = Am + n is 
actually the most general nontrivial one, dictated by the linear A-dependence in (22). 
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The major gain here hes in the fact that there is a systematic expansion of the phase in powers 
of 1/A, while the leading term is of the order 0(1). As it will turn out, this will become very 
important as we go on to study the asymptotic behavior of the whole wavefunction. We shall 
systematically calculate everything up to terms of the order 0(1/A^), since this is the lowest 
self-consistent approximation for the Laplace method, as we shall see below. This is also the 
reason why we keep the 0(1/A) term in (22). 

5 Stationary-point approximation 

We have written the wavefunction (28) in the form required for the application of the Laplace 
method. However, given that the integral is multidimensional, there are certain complications. 
The first step is to expand the phase into a power series around an extremal point y^. We have 
to cast the series in the form which separates the variables y^, which are not to be integrated 
over, from the internal i/n {N G I) variables 

S{A, y) = S{yl, y*) + J^ 7^- (l/iv - y^) + X] 7^ (^i " ^i) 



N£X 



dyN V ^Vj 



4E 



ivM - yh) ivN -y*N)+ Y^ 75 — jf- ivM - yh) {yj - y]) 



2 M%x ^VMdyN ^^ dyudy, 

3 



-IE 



{vj - Vj) {yr - y]') + ■■■ 



2 j-f dyjdyf 

Given that all derivatives above are evaluated at an extremal point y^, the terms with first 
derivatives vanish. Also, since all the differences y^ — y*^ go to zero as 1/A, due to equations 
(26) and (27), we can neglect the terms of the order {y — y*)^ and higher. This leaves us with 

SiA, y) = S* + ^iy- y*fS"iy - y*) + (%■ - y*fS\y - y*) + ^{y, - y*fSiy, - y*), 
where we have introduced a shorter matrix notation. 



S* = S{A,y*), S" 



dyMdy_ 



N 



s' 



d'^s 



s 



d's 



dyjdyu y*' dySy 



j^yy 



{M,N el). 



The matrices S" , S' and S are of the type uj x uj, J x uj and J x J, respectively. At this point 
we see that keeping all terms of the order up to 0(1/A^) is necessary, since if we had kept only 
terms up to 0(1/A^), the phase would have been approximated by a constant, and the Laplace 
method would have not worked. 

Now the integrals in the wavefunction obtain the form 

d^y eA'^(A'J^) 

^ gA2S*g^(%-?/*)^5(j/,-yj) f ^^y ^^(y-y')Ts"iy-y*)+K^(y^-y*)Ts'{y~y*) ^ 

At this step it is important to note that we have one integral of this type for every allowed value 
of m^ and n^, and for every extremal point y*^ which is inside the integration domain D'^. The 
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Gaussian form of all the terms in (24) guarantees that all extremal points of the integrand are 
maxima, which in turn means that all eigenvalues of the matrix S" are negative or zero. 

This allows us to expand the integration domain D'_^ to M.^, since the eventual "exterior" 
extremal points are not taken into account while everything else is negligible in the limit A — )■ oo. 
The zero eigenvalues contribute with linearly divergent terms, but this can be regularized in the 
sense of the generalized Gaussian integral (69) (see Appendix E). Even when the integration 
domain is extended to M^ the integral converges (or has a constant divergent contribution), and 
the integration can be explicitly performed by using the formula (69). 

Therefore, after a suitable orthogonal change of variables z = 0{y — y*) which brings S" and 
S' in a block-diagonal form, we perform the integration and obtain 



y*€D' 



'si 



yj-Vj 



dz 




-(yj~y*VNM-^N'^{y,-y*) 



Here r is the rank of matrix S" , while M = Mg" (as defined in Appendix E). It is important to 
note that the orthogonal change of basis which brings S" into a block-diagonal form does not 
necessarily guarantee that S' will also reduce to zero in the null-space of S" . In other words, if 
K denotes the null-space projector of S" , it is not guaranteed that S'K = 0, which was assumed 
in the above equation. However, if this assumption is violated, we can immediately conclude 
that ^_(_ does not have the desired asymptotics, because of the integral of a linear exponential 
function. Nevertheless, numerical investigations (see subsections 7.1 and 7.2) suggest that S'K 
is indeed zero, so we proceed assuming that this is satisfied. 



Using the fact that A{yj — y* 
write the resulting integral as 



j* we can switch back to the variables j and jo, and 



E 



^A^S* 



A^ 



(2vr) 



where we have introduced the matrix 

S = 



detM 



dz 



-(x-j*y six-r) 



s - nm-^n'^ 



(31) 



This matrix is of type J x J, and represents the key ingredient of the calculation. It is known as 
the Schur complement and is well-studied in general matrix theory (see Ref. [11] and Appendix 

F). 

Finally, we substitute this result back in the equation (28), and obtain 



*+(j,io) 



2^A^+^-^ 



\/l2^i?4(0,e- 



V 



(32) 



E 



mi,ni<CA 



■ E 

mv^nvGl' y*&D' 



y^ ('_l)Et,(Ami,+n^) 



„A2 5* 




dz 



'ij-j*y sij~j*) 



The matrices M, N and their rank r typically depend on the initial choices of parameters 
TTiv, Uv and y^ G D'_^_. Due to the nature of the Laplace method, we should keep in the sum only 
those parameters which give minimum r, and among those only the ones which give maximum 
S* . The resulting sum of the remaining Gaussian functions represents the asymptotic behavior 

of *_L. 
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6 Asymptotic behavior of the matrix S 



Let us now return to (28), 



^+(j,io) 



A \A+Lu-- 



2^A 



\/l2^i?4(0,e~'^ 



E 



mi,ni<CA mv,n.vCA 



Jd' 



The phase 5'(A, y) can be expanded into a power series around an extremal point y^ (the first 
derivatives vanish, while the third and higher-order derivatives are of 0(1/A'^)), so that 



SiA, y) = 5(A, y*) + -{y - y*)' A(y - y*) 



Here the A matrix is defined as 



A 



MN 



d^S{K,y) 



dyudyN 



y=y* 



It is of type $7 x fi, and it is convenient because it can be decomposed into blocks of size J and 

CO 



A= [A 



MN\ 



s 


S' 


{sr 


S" 



After an orthogonal transformation of the basis, the matrices S" and S' will reduce simul- 
taneously into a block-diagonal form, so that the A matrix will obtain the following form 



A 



s 


N ' 





M 
. 



(33) 



Integration over the zeroes in (33) will boil down to a trivial divergent part, as seen in (32), and 
we can consider only the nonzero block. Note that S matrix defined by (31) is actually (minus 
one half of) the Schur complement of the nonzero block of the A matrix. 

At this point we apply the theorem from Appendix F, which states the following: 

• R = r + p, where R, r and p are ranks of matrices A, S" and S respectively, 

• if p = J then det Ma = ± det M det 25*, 

• if < p < J then det MA(det B^f = ± det M det M^g, 

where the signs depend on even/odd rank of S. See Appendix F for a proof of the theorem and 
the definition of the matrix i?4. 

The A matrix can be expanded into power series 

A = Ao + ]-Ai + ^A2 + ... , 

and as we shall see in the next section, the first nonzero leading term in the series is always Aq. 
Consequently the determinant for the A matrix is given by 



det A = det Ao + O 
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Given that M is a subniatrix of A, it follows that 

det M = det Mo + O 
Assume now that 



where n G No- If the rank p of S* is positive, we have 

M25) = ^det(25„) + 0(^). (34) 

Equation (34) implies three distinct possibilities. If p = J, then we can use the first identity 
for determinants from the theorem, and obtain the equation 

det Ao = i-r;^ det M det(2S'„) + o(- 

which is consistent if and only if n = 0. If < p < J, we can use the second identity for 
determinants, and obtain the equation 

detAo(detfi4)^ = ±— ^ det Mdet(2S'„) + O [- 

which is consistent if and only if n = 0, due to the fact that deti?4 ^ 0(1) (see remark 3 in 
Appendix F). Finally, if p = we have S = 0. 

Therefore we have essentially two possible situations. If p = 0, then the matrix S is equal 
to zero, because p is its rank. This implies that the wavefunction (32) is constant in the leading 
order of A, i.e. the j dependence only appears in the subleading terms, which are of the type 
(j ~ io)^ a-iid higher. On the other hand, if n = 0, 5 is different from zero in the leading order, 
so that (32) is a Gaussian, but not of the required type (4). Hence the large-spin asymptotics 
of '1'+ is never a Gaussian function of type (4). 

7 Computation of the matrix S 

The main result of the previous section has been obtained under the assumption that the matrix 
A has a leading contribution of 0(1) in A. In this section we will demonstrate this by an explicit 
computation for some concrete spin-network diagrams. In order to do so, we need to explicitly 
find an extremal point y^ and evaluate the corresponding A for a given spin network. 
An extremal point is defined by 

0, wherey*eA', NgM, (35) 



dvN 



where we must also take into account the consistency conditions (26) and (27). 
Differentiating (29) and (30) we obtain 

a5„ i|5, i|s,^ /ly (36) 



dyN ^dvN ^^ dvN \^-^ 
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dSo 
dvN 
dSi 
dvN 

dS2 
dyN 



-2C{yN - yo)SN,x 



y^ On,v ^ ys0s,v - w-„7r , 

Ndv 



Yl ^^.- (i 



n„7r 



VN V [2K, dyN \i^v I 



dFpL 
dyN 



where we have used the Schlafli differential identity (56) for a tetrahedron 






Vs 



dy 



0, ViV, V . 



N 



We will also need second derivatives: 
dyudyN 



dyudyN 

d^S2 
dyudyN 



-2C6mnSn,\ — / ^ 

V 

M,N&v 

Ed9N,v /TT 
dyM V4 

Vn V 

M,N£v 



Om,vOn,v + — — - I 22'^"^"''" ~ ^""^ 



dyjvi 



\S(^V 



1 92 K 



+ — rDL + 



(37) 



+On,v^ 1- (^M,v 



2K, dyudyN 2V^ dyM dyN dyM 

d^FoL 



\s6f / 



a^M ayAT yj^ J dyMdyN 

A solution of the system (36) can be obtained perturbatively in 1/A via the ansatz (25), 

Bn , Cn „ / 1 



y^ = A^ + — + ^ + O 



In the lowest order we obtain a nonlinear system for A 



N, 



2C{An - yo)SN,x + Yl ^N,v{A) ( Yl ^sOsA^) -m^A = 0, 

V \s£v J 



which can be reduced to a simple equation using (26), 

^A = ya- 



rn 



This was already guessed before based on the analysis that the extremal point must be in the 
vicinity of extremal points of the cosine and the Gaussian functions ^{x\). 

At 0(1/A) order we obtain also a linear system of equations for the Bn coefficients 



E 

M 



2C5mn5n,x+ V 5^1 ( VAA,.(v4)-m,7r| + V 0jv,.(^)^M,.(A) 

M,Nev M,Nev 



B 



M 



- Y ^nAA) 



V 

Nev 



TT 



n„7r 
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which can be siniphfied using (26) and (27) to 



(39) 



Similarly, by applying (26) and (27) at the 0(1/A^) order we obtain a linear system for the 
coefficients Cn 



E 

M 



'^C6mnSn,x + 2J (^N,viA)dM,viA) 



M.N^v 



Cm 



(40) 



An V K^VvdyN a ^ ' ) 



Let us now write the complete set of equations which determine the extremal point y* 
• equations for A^: 



y^ Ases^,a{A) = m^TT, Ax = 2/0, 



sSf 



• equations for B^: 



equations for Cn'- 



^Bses,v{A) = n^7r-j, Bx = 0, 



(41) 



(42) 



s£v 



M 



2C6mn6n,x+ Y1 GN,v{A)eM,viA) 

V 

M,N&v 



Cm 



(43) 



uM--T. 



An 



1 dV„ 



-;f \2VvdyN 



+ eNAA)FDLiA) 



Before we engage in finding solutions to these equations, let us introduce some notation and 
discuss the form of the second derivatives of the phase S{A,y). We will introduce the following 
shorter notation for various derivatives evaluated at the particular extremal point: 



On,v = On,v{A), 9mn, 



dON,v 



d'^dMv 



dyM 

K = K(^), Vn, 



TMNs,v 



A' '" dysdyM 



SMNpn,v 



d^e 



N,v 



dypdyqdyM 



dyN 



Vmn,i 



d^V, 



dyudyN 



By using (37) and the ansatz (25) for y* , as well as the equations (26) and (27), we obtain 
f 



A 



MN 



2C6mnSn,x + 2^ (^m,vOn,v 



V 

M,N<^v 



1 

A 



7 ^ 7 ^ Bs {9m,v(^Ns,v + 9n,v0Ms,v) 



^M,N£v 



1 

'A2 



--T2~^MNSN,a + 2^ 



^N 



V 

M,N£v 



+ OmN,vJ^DL + ^N,: 



2K, 



2K? 



dyM 



, . dFoL 

+ ^M,v^^ 

A dyN 
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+ E 

V 

M,N<^v 



Om,v > , CsOns,v + 0N,v 2_^ CsOms,v + ^^M,?) 2_^ BpEgO^pg^y + -9n,v 2_^ BpBqOM 



set) s£t) PtQ&v p,g&^ 



+ 2^ ^MTV.i 



/ ^ C's^s,?; + ( 7 ~ '^t'Tl" j 2^ Bs9mNs,v + - 2^ BpBqOpg^y I "*" ^ ( Xs ) 
,sGt> sSn Pi9G^ J 



• (44) 



M,Afei> 

Therefore, after finding an explicit extremal point, we have substituted the coefficients A]\[, 
B]\[ and Cjv into the above equations and obtained an explicit expression for the A matrix. 
The curly braces group all terms of orders 0(1), 0(1/ A) and 0(1/A^), respectively. This 
demonstrates that the leading order of the A matrix is an 0(1) term. 

Given that all of the above equations are fairly complicated, we will investigate them order 
by order in A. 

7.1 The 0(1) approximation 

Let us rewrite the equation (41) in this approximation as 

Y^ AsOsA^) = m„7r, Ax = yo. (45) 

We will also rewrite the expression for the A matrix as 

Amn = -2C5mn5n,x- Y1 ^mA^WnA^)- (46) 

V 

M,Nt^v 

The calculation of the matrix S can be organized in the following way: 

• calculate all extremal points y^ = A^ by solving equation (45); 

• for every extremal point obtained, calculate A matrix according to (46); 

• split the A matrix to blocks S, S' and S"; 

• determine the rank of S" , and the null-space projector K if the rank is less than w; 

• check whether or not S'K = 0; if it is nonzero the procedure fails and the wavefunction 
does not have Gaussian form; 

• determine an orthogonal matrix O which diagonalizes K and use it to change to a basis 
where matrices S" and S' are block-diagonal; read-off the nonzero blocks M and TV; 

• compute the matrix S according to the equation 

S = --(S-NM-^N^ 
2 V 

The system of equations (45) is highly nonlinear, and therefore extremely hard to solve. 
Nevertheless, one exact solution can be guessed 

yj^ = An = yQ, On,v = Oq ^ tt - arcsm —- = arccos ( - o ) ' "^^ = Vo- 

We will refer to this solution as the "diagonal" solution. It corresponds to a situation 
where all the tetrahedra are equilateral, which is a very symmetrical configuration. This greatly 
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simplifies the equations, and the required matrices can be calculated without additional approxi- 
mations. In the case of the loop spin network and the theta spin network a computer calculation 
gives 

S'K = 0, R = r, S = 0. (47) 

This result agrees with the p = case discussed in section 6. 

The calculation of the S matrix can be also performed numerically for non-diagonal solutions. 
In fact, in the numerical approach, the hardest first step, solving the equation (45), can be 
completely sidestepped. This is due to the fact that A depends on Ai\f only through the angles 
6]\f^y{A). Thus one can design an algorithm which chooses the angles completely randomly from 
their domain [0, vr], which in principle covers also the angles obtained using any specific solution 
of (45). The constant C > can be also chosen randomly. 

One such algorithm has been implemented on a computer to explicitly calculate the A matrix, 
the projector K, the rank of S" , then S'K, M and N, and finally S. It has been executed 50 
times with random initial data for the cases of a loop spin network and a theta spin network, 
and each execution gave the same result (47). The numerical precision of the calculation was 
10^^. The precision can be arbitrarily increased at the expense of the execution time, and in 5 
executions the precision was raised to 10~^^, with no change in the result. 

7.2 The 0(1/A) approximation 

Let us rewrite the all necessary equations up to 0(1/A^). The extremal points can be found 
using the ansatz 

y^ = Aat -F — . 
The coefficients A^' and B]\f will be determined by the equations 



Y^ As9s,v{^) = rriyTT, Ax = yo , (48) 

Y^ B,es,M) =n,TT-j, Bx = 0. (49) 



> , (50) 





sdv 




A matrix is given as: 


r 


^MN = -'^C5mn^n,\ - 


- 2^ 0M,v0N,v - -T ' 


1^ l^BAt 




M,N(^v 


_M,Nev 


where 


r\/-i 




On,v = dN,viA), 


a O0n,v 

oy. 



The procedure for computation of S matrix is the same as in the 0(1) approximation, up to 
two additional steps. These two steps consist of solving the linear system of equations (49) for 
-Bjv coefficients, and then expanding the resulting matrix into a power series in 1/A, 



5 = So + f + o(^) 



This case can be also analyzed analytically (by using the diagonal solution), as well as 
numerically (by randomly generating the A^r coefficients). In contrast to the previous case, one 
should also always specify the n^ parameters as a part of the initial data, in order to solve (49). 
One algorithm for this has also been implemented, and we obtained the same result (47) for a 
loop and a theta spin network. 
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8 ^_ integral 

The asymptotic analysis of ^_ can be done by using the same method as in the ^+ case. 
However, it turns out that ^_ analysis is considerably simpler due to qualitatively different 
nature of the scaling laws for the asymptotic expressions for 6j symbols when some of the spins 
are not large. 

The ^_ part of the wavefunction can be approximated by an analogous expression to (20) 

z -'- a a A t) 

where the vector z denotes the small spins, while the vector Ay denotes the large ones. We 
will sum over the small spins rather than integrate and the number uj of large spins is smaller 
than in the ^-|- case, while F(z) represents the part of the wavefunction which does not depend 
on large spins Ay. Which spins can remain small and which must be large depends on the 
triangle inequalities built into the 6j symbols. These restrictions on the spins z and Ay will also 
depend on the detailed topology of the spin network. However, these details will not affect the 
asymptotic analysis. 

The next step is to use the asymptotic formulas for the 6j symbols, but now we cannot use 
only the PR formula (57) since the 6j symbols where not all of the spins are large will appear. 
The asymptotics of 6j symbols with 3, 4 or 5 large spins is given in the Appendix, see formulas 
(59), (60) and (61). 

Let the asymptotic behavior of each 6j symbol be described by a function (j)i{Ay, z), so that 

{6j{Ay, z)} ^ (t>i{Ay, z) as A ^ oo , 

where i = 3, 4, 5, 6 denotes the number of large spins. When i < 6, we will refere to the 
corresponding 6j symbols as degenerate. When i = 6, we will refere to the corresponding 6j 
symbol as a non-degenerate. 

Note that in a degenerate case 

cj)i{Ay,z)^ AP' fi{y,z) as A ^oo, (52) 

while in the non-degenarte case 

(/>(Ay) ^ AP cos (A/(y)) as A ^ cx) . (53) 

Let us write the integrand of ^_ integral in an exponential form suitable for stationary point 
approximation. The corresponding A matrix is given as before 

_ d''S-{Ay,z) 

oyMoyN y=y* 

where the new phase S- is now proportional to ^^ ln(/)j,(Ay, z). It is easy to see that the A~ 
matrix will contain the terms of the form 



1 



^v 



dAyJ 

and 

1 5V. 



(pv dAydAy 

All these terms are of 0(1/A^), since the asymptotic functions (pv behave well when dif- 
ferentiated by Ay, see (52). Furthermore, given that the A matrix should be evaluated at an 
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extremal point y*, one can see that the 6j symbols which contain only 3 or 4 large spins give 
a sub-leading contribution, since they decrease monotonically in the limit A — )■ cx) and thus do 
not have any extremal points. The 6j symbols with 5 large spins have extremal points, but they 
also provide only a contribution of 0{1/ h?) to the S^ matrix. The 0(1) contribution to S- 
comes only from those Qj symbols where all spins are large, due to the fact that the argument 
of the cosine function in (53) is proportional to A. 

Therefore the ^_ integral will have the same asymptotics as the ^-|_ integral if there is at 
least one non-degenerate 6j symbol contributing to S- . If only the degenerate 6j symbols are 
present, the contribution to S'_ will be of 0(1/A^). 

9 Conclusions 

We have shown that the leading term in the large-spin asymptotics of a flat-space wavefunction 
^fc(r,io) = ^(i, jo) is given by 

*(j,Jo)»2^iV„p(jo)e ^0 +2^N+g{jo)e20jo)B„,ujo)^ (54) 

n,p n,q 

where n £ Z, p £ S-, q & 5"+ and B and D are constant (independent of jo) matrices. The 
index sets S± correspond to the stationary points of the large-spin approximations f± of / in 
the regions D± and the corresponding functions A^ (jo) will be powers jg*, r± G Q. In the case 
of a loop and a theta spin network the computer results are consistent with B^ = C^„ = for 
all n and all q which means that in those cases the matrix S vanishes. 

The asymptotics (54) is not of the type (4) required for the correct semi-classical limit. This 
means that the wavefunction ^fc(r) jo) will not give the correct graviton propagator asymptotics. 
One can argue that some other wavefunction may give the correct asymptotics, but the problem 
is to see what other wavefunction can replace ^fc(r, jo). Although our result applies only to 
the Euclidean LQG, one wonders what is the relevance of this result for the Lorentzian LQG, 
given that there is a strong belief that the Euclidean and the Lorentzian theories should be 
related by some kind of an analytic continuation. Note that Lorentzian analogs of the Euclidean 
wavefunctions used in this paper are not known. However, one can try use one of the recently 
proposed Lorentzian spin foam models [12, 13, 14] in order to construct a Lorentzian spin- 
network wavefunction. The large-spin asymptotics could be then computed by using essentially 
the same techniques as the ones introduced in this paper. 

As far as the our result is concerned, there are certain caveats. The first caveat is that the ob- 
tained asymptotics is for the Ponzano-Regge regularization ^^ (F, jo) of the zero-cosmological 
constant spin-network wavefunction ^(F, jo). The wavefunction ^|j'^(F, jo) is different from the 
quantum group regularization ^^(F, jo) of ^(F, jo), but the physics intuition suggests that the 
asymptotic behavior of ^^ and ^^^ for large k should be essentially the same, up to constant 
factors, see for example [15]. However, it still remains to be proven that these two wavefunctions 
have the same semiclassical asymptotics. 

The second caveat is that the construction of ^^ wavefunctions is triangulation dependent 
and our calculations have been done for the simplest triangulation of S^. However, it is not 
difficult to see that by taking a more complex triangulation, the corresponding asymptotics 
will not change qualitatively, because the corresponding f{x,j,jo) will be always a product of 
6j symbols and our method of computing the asymptotics is independent of the number of 6j 
symbols. Similarly, one can hope that a special choice of the function ^(A) can lead to a desired 
asymptotics. However, a generic /i(A) gives a contribution to S of 0(1) in jq. Note that it is 
possible to fine-tune fi for a given spin network such that the 0(1) contributions cancel and the 
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O(l/jo) contributions are non-zero. However, this fine-tunning depends on the spin network 
and hence one cannot find a /x which will work for all spin networks. 
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Appendix 

A The spin-network diagram 

The spin- network whose evaluation appears in (11), is given in the case of a single loop spin 
network with spin j by the following diagram 
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Figure 3. 
All spins take values from the set {0, |, 1, |, 2 . . . , 2}- Here n is the degree of the corresponding 



,iTr/k+2 



quantum group SU{2)q such that gr^^+^ = 1 i.e. q = e* 

The evaluation of the above diagram, to which we refer as the amplitude, can be calculated 
via the following rules: 
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The amplitude of the whole diagram is the product of the amplitudes for the vertices. 

The amplitude of a three-vertex is proportional to the amplitude of the corresponding 
tetrahedron spin network, as follows 




= M 




The amplitude for a four-vertex is proportional to the amplitude of the corresponding 
tetrahedron spin network, which is given by 




■M 




The normalization A^ is given by, see [4] 

N 



v/|e(a, h, c)Q{c, d, e)e(a, e, /)e(6, d, f)\ 

where ©(a, b, c) is the evaluation of the 6'-graph 

n( h ^ - ! -,^a+h+c [a + b-c]q\[a + c-b]q\[b + c-a]q\[a + h + c+l\q\ 
""^"'"''^-^ '^ [2a],![26],![2c],! 

B Tetrahedron spin network and the 6j symbol 

The amplitude for a tetrahedron spin network is, roughly speaking, the value of the correspond- 
ing 6j symbol. Following the conventions of [3] and [16], the amplitude is given by the equation 




a h c 
d e f 



Vl©(a, b, c)eic, d, e)e(a, e, /)e(6, d, f)\ . 



(55) 
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The tetrahedron associated to a 6j symbol has the fohowing geometric properties. The 
length of an edge colored by a spin j, is given by 

•^ 2 
The area of a face, with edge lengths li,l2,hi is given by the Heron formula 



A 



123 



^/s{s-li){s- 12){S - h), 



where « = -(^i + ^2 + ^3) • 



The volume of the tetrahedron is given by the Tartaglia determinant 
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where lij is the edge length between the vertices i and j. Note that if we fix all 6 edge lengths, 
we could create in total 6! different tetrahedra. However, some of them are equivalent up to 
rotations and reflections. When this is taken into account (4! permutations of the 4 vertices of 
a tetrahedron), we end up with the total of 6!/4! = 30 possible inequivalent tetrahedra. Those 
30 tetrahedra have different volumes, which means that if we provide six numbers lij, there is 
30 inequivalent different ways to position them in the determinant above. 
A dihedral angle of a tetrahedron is given by the formula 



sin( 



3 aVgbcdef 
2 ^abf^ace 



The angle 6a corresponds to the edge a, and it is constructed between the outward-normal 
vectors of the faces a, b, f and a, c, e, such that the edge d does not intersect the edge a. This 
angle is equal to tt — (fa, where ipa is the angle between the faces a, 6, / and a, c, e. 

Regarding the above equation, once we fix the value of the right-hand side, there are in general 
two different angles which satisfy the equation. This is a consequence of double-valuedness of 
the arcsine function on the [0, tt] codomain. However, when the edge lengths are given, 9a can 
be expressed as 



TT 



arcsm 



3 aVr. 



abcdef 



2AabfA, 

'3 aVn 
arcsm I — ■ 
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abf^ace 



if d^ 
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Finally, it is convenient to notice that the angle 9a is unchanged if we scale all edge lengths 
of the tetrahedron, i.e. if we multiply them all with the same positive constant. This is easily 
visualized since the tetrahedron does not change its "shape" if we "zoom in/out" , and can be 
verified analytically by inspecting the above equation. In other words, we have an identity 



9Aa(Aa, Xb, Ac, Ad, Ae, A/) = 9aia, b, c, d, e, /), 



A>0. 



By differentiating this identity with respect to A and a, b, c, d, e, f, it is easy to derive the 
so-called Schldfl,i differential identity for a tetrahedron 



Y-7 ^ 
^, dl, 

s=l '^ 



0, yipG {a,b,c,d,e,f} . 



(56) 
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As far as the 6j symbols are concerned, we are mainly interested in the asymptotic formulae 
for large spins. The limit of large spins is defined as a limit when the scaling parameter A tends 
to infinity. The relation between a spin and A is defined by 



hm fij) ^ hm f{Ak - I: 



where 



h 



..1. 



Here k is some initial finite spin, which is scaled to j via A according to the map 

j + - = a( k + - 

In other words, it is not spins themselves that are being scaled, but rather their corresponding 
edge lengths. 

The large-spin asymptotics of a 6j symbol can be defined in the following way. Fix the 6 
initial spins ki, . . . ,ke, and associate to them a tetrahedron with the edge lengths k = ki + ^. 
This tetrahedron is then scaled by the parameter A into a new tetrahedron with the edge lengths 
AZj. Let us denote the corresponding spins as ji, then the following asymptotic equation holds 
[9, 10]: 



lim <! ^' -^.2 -^3 ;> ^ lim 

A^oo 



j-^00 [ J4 J5 JQ 

where 



Ah - 2 
A/4 



1 Al 
2 



1 

2- 2 

5- 2 



Ah 



1 A/,-1 A/fi 



y/l2TTA^V{l) 



cos5(A, /) , 



5(A, /) = AJ2lses{l) + J + jFnLil) + O (j^) 



(57) 
(58) 



The first term in S represents the familiar Regge action, while F£,l{1) is a very complicated 
correction of 0(1/A), see Ref. [10], equation (32). 

The second possible configuration is when 5 spins in a 6j symbol are large, of 0(A), while 
one is of 0(1). In this case, we have the following asymptotic formula 



lim 



Ji 



J2 



J3 



{-l)n+n+T^+k2+k, 



dZL iO) , 



where 



j^oo ( J2 + k2 ji + ki ks j \/{2jx + I)(2j2 + 1) ''''''' 

ii(ii + 1) + i2(i2 + 1) - is (is + 1 



(59) 



cost 



2Vii(ii + i)i2(i2 + i) 



^ 61 ^ TT. 



Here d^ ^^ iff) are the usual matrix elements of the SU{2) rotation operator, see [9]. 

Next we have a configuration with only 4 large spins, and the corresponding asymptotics is 
given by 



lim 



h j + k2 j + /C3 

ki j + ^5 i 



_'y\ki+k4+mm(k2+k5,ks) 



(2j)- 



-1-|A;2+A;5-A:3| 



\k2 + k5 - /csj! 
{ki -k2 + k^)\{ki - k^)\{k^ -k5 + ^3)1(^4 - k2y. 



1+0(4 



isgn{fc3-fc2-fc5) 



(60) 



_iki +k2- ksV-iki + /C5)!(/C4 + fes - ksYiki + k2)\_ 

see [9]. 

Finally, the configuration with only 3 large spins has the asymptotic formula proportional 
to the Wigner 3i symbol [9] 



lim 



ki k2 

i + A;4 i + ^5 



^3 1 ^ / 2)^'l+fc2 + fc3+2(fc4+fc5) ^ [ 

i / ^/27V 



ki k2 k-i 

/C5 /v4 i\'4l "'S 



(61) 
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As a final remark, note that the spins and the corresponding edge lengths for a 6j symbol 
are constrained by the triangle inequalities. However, these constraints are not sufficient to 
guarantee the condition V^ > for the volume of the 6j symbol tetrahedron. There are choices 
of spins such that the 6j symbol is well defined, and the left-hand side of (57) is real, while 
V'^ < 0. Consequently the right-hand side of (57) is a complex number. Of course, in such 
situations the equation (57) does not apply, and we have a different, exponentially decreasing 
asymptotics, see [9]. However, we are interested in the extremal points of the asymptotic formula, 
and the exponentially decreasing asymptotics does not have any such points, so it is not relevant 
for our purposes. 



C The cosine approximation 

In this section we shall give a rigorous derivation of the "cosine formula" (23). We start by 
introducing the Heaviside step-function as 



H{x) 



for X < 0, 
for X = 0, 
for X > 0. 



This is a standard definition, and the standard rules apply. For example, derivative of a 
Heaviside function is the Dirac 6 function and the integral of 5- function from — oo to x gives 
H{x). We can then define the so-called rectangle function 



Kaix) = H{x + a)H{-x + a) 




1 

2 
1 
1 

2 





for X < —a , 
for X = —a , 
for — a < X < a, 
for X = a , 
for X > a . 



The function Ka{x) is equal to one in the interval (- 



-a, a) 



zero outside, and at the points 



-a and a K is conveniently defined to be ij, so that it has a nice property proven in Lemma 1. 
Lemma 1. The following identities hold 



Kk{x + (a + b)) + Ka{x) + Ki,{x - (a + b)) = K, 



'a+2b{x) , 



Y^Ka{x-2ka) = 1. 



(62) 
(63) 



kez 



Proof. The first identity can be demonstrated by using the definition of K. First, we see that 
outside of the interval [—a — 2b,a + 26] both the left-hand and the right-hand sides of (62) are 
equal to zero. Next, the three terms on the left-hand side of (62) are equal to one respectively in 
the intervals (—a — 26, —a), {—a, a) and (a, a + 26). Finally, at points —a and a the appropriate 
terms are equal to 2 ^ind add up to one so that the resulting function is continuous and equal 
to one in the whole interval {—a — 26, a + 26). At the boundary points of [—a — 26, a + 26], 
the left-hand side of (62) gives the contribution of ^. All these results taken together form by 
definition the right-hand side of (62). 

The identity (63) can be proved by applying the identity (62) iteratively. First notice that 
the sum over all integers is in fact defined as the limit rri — )■ 00 of the sum over the domain 
—m, —m + 1, . . . ,m — l,m. Consequently we can write: 



fcez 



Ka{x-2ka) 



m 
lim y 

k=—m 



K„ 



2ka). 
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Now note that by applying the first identity m times we have 

m 
^ Ka{x- 2ka) = K(2m+l)a{x) ■ 
k=—m 

This identity can be also easily seen graphically, and can be proved by induction over ra. There- 
fore, 

y^ Ka{x- 2ka) = lim i^(2m+i)a(a;) = Koo{x) = 1 . 

fcez 

End of proof. 

The function Ka{x) was introduced because it allows for a neat "cutting" of the appropriate 
pieces of the real line. Now we make use of this property in order to prove the following identity. 

Lemma 2. The identity 



cosx = N^(— 1) Kil{x — kT:)cos{x — ki:) 



kez 
holds. 

Proof. Start from the right-hand side and compute the left-hand side in the following way 

RHS = > (— 1) i^ 2L (x — fcvr) (cos a: cos fc-TT + sin X sin /cvr) 

fcez 
= ^(-l)'=K|(x-fc7r)(-l)'^cosx 

fcez 
= cos X y^ Ke. (x — kn) 



2 

fcez 



cosx . 



End of proof. 

Now if we note that Kn^x — kn) cos(x — kir) is continuous and non-negative for all x G M and 
k £ Z, the statement of Lemma 2 can be rewritten in the form 



cosx = 

fce 



Sr^ / _'^\k ^^ii[K :^{x-kTT) cos(x-kTr)^ ,g^^ 



which represents the key "cosine formula". The exponent can be expanded into a power series 
around the point x = /cvr, 



i^i(x — /cvr) cos(x — /cvr) = — (x — fcvr)^ [x — kn)^ + 0{x — kir)^, (x — )■ /cvr) , 



2 



In 
and in the leading order approximation we can write 

COSX« J](-l)V5(--'=-)'+^(-''=), 



kez 



where R{x, k) is the remainder of 0((x — /ctt)^) when x — )■ fcvr. 

This kind of approximation is useful since the cosine function is well approximated in the 
vicinity of all extremal points simultaneously. In the crudest approximation, the remainder 
i?(x, h) can be substituted by some average value over one period of the cosine, R, so we write: 



COSX«e-^J^(-l)V5(^~'=-)', X^kTT. 



fcez 
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The constant R can be calculated, because at any point x = fcovr (/cq G Z) the formula must 
become exact, with no error. Namely, expanding the exact exponent into a power series and 
discarding certain terms we effectively allow for the appearance of the "tails" in each Gaussian 
in the sum. These tails then give artificial contribution to other Gaussians, and the constant 
e accounts for the appropriate "correction" . For example, substituting x = into the left and 
right side, we have 

kez 

so that 

e^ = — « 1, 01459 . 

^4(0, e-V) 

Here 'd4{u,q) is the so-called inverse elliptic theta-function of the fourth kind. Of course, if 
we had kept the x'* term in the exponent, the Gaussian tails would have been different and the 
constant R would have had difi^erent numerical value, but again such that the equation is exact 
at all points x = kn. If we had kept all terms in the exponent, the tails would have vanished 
and this constant would have been equal to one. 

Anyway, we are interested only in the crudest approximation of the cosine function with 
Gaussian functions, so that the formula we will use is 

cosx« L__^(_i)fce-5(^-M', x^kir. (65) 

^?4(0,e~V) ^g2 

One very important remark here is that this equation makes sense only in the vicinity of 
cosine extremal points, i.e. in the neighborhood of the points x = kir. Therefore one should make 
sure that the argument of the cosine is close enough to some integer multiple of vr, whenever one 
applies the formula. For the points in the vicinity of kiT/2, the approximation is less accurate. 

D Laplace's method 

Let f{x) be a real function in the interval [a,b] such that it has a single global maximum at a 
point xo G [a, b]. We would like to find the asymptotics of the integral 

1= I dx e^^(") , 



for large A. We will describe here the Laplace method of calculating this asymptotics. 

The idea of the method is the fact that e^^^> goes to infinity at the fastest rate at the point 
of the global maximum as A — )■ 00. One can then approximate e ■'^^' in [a, b] by a Gaussian 
function centered at xq. This also implies that the integration domain [a, b] can be extended to 
the whole real line, since the "tails" of a Gaussian function do not give a significant contribution. 
By using the Taylor series for /(x) at the point xq we obtain 

/(^) = E ^^^-^(^ - ^0)" = /(^o) - l\f"i^o)\ix - xof + 0(x - xo)^ . 

n=0 

By changing the integration variable to y = x — xq, the integral / is reduced to a sum of 
Gaussian integrals 

-' 00 ^„ \ 



./hi \ — i n\ 



\n=0 



27 



which imphes 



=,A/(a;o) , 



27r 



l + O 



(66) 



V Mf"{^o)\ 

In the case when there are several global maxima in [a, b] with the same value of f{x), then 
the formula (66) becomes 



/ 



xo£[a,b] 



pA/(xo) 



2tt 



A|/"(xo)| 



1 + 



(67) 



E Generalized Gaussian integrals 

Let us start from the identity 



^3,g-|aa;2+6x 



e2« , a > . 



This identity can be generalized to 



d"xe' 



\x'^Ax+Bx 



(27r)" 

det^ 



JBABT 



(68) 



where A and B are matrices of the type nxn and 1 x n respectively and A is a symmetric matrix. 
This identity can be derived by reducing the integral in (68) to a product of one-dimensional 
integrals by a change of variables x' = Ox, where O is the orthogonal matrix which diagonalizes 
the matrix A. In the diagonal basis the integral reduces to a product of n one-dimensional 
integrals. The critical assumption here is that the eigenvalues of A must be strictly positive 
if the integral is to converge. However, we are interested in the situation when some of the 
eigenvalues of A are zero. In this case the integral in (68) diverges, and our goal now is to 
regularize this divergence in one special case. 

Let us start with an integral over a compact domain D = [— G, G]" 



Id 



cTxe 



-^x'^Ax+Bx 



D 



For simplicity, assume that the matrix A is already diagonal, and denote its rank as r < n. Also, 
denote all n eigenvalues as a^, and let ai, . . . ,ar / 0, while Or+i, • • • , a„ = 0. Next, let K be 
the projector to the null-space of matrix A, and assume that BK = 0. This means that in this 
particular basis we have Br+i, ■ ■ ■ , Bn = 0. Now the integral //) can be split into a product of r 
Gaussian one-dimensional integrals and n — r one-dimensional integrals of a constant function, 
where the constant is equal to one, 

r-G , ^ f-G 



fo=n 



dxi 



' O.'^i ^ii-'-^i^i 



-G 



n 

;=r+i 



dxj 



2 "^i ^^ i-^i-^i 



-G 



In the limit G — )■ oo we will define a regularized / 



hmc^oo Id 



lim 



G^oo 



rG 



f_^dx 



(2vr)'' ijsfM-^NT 



detM 



(69) 



Here M and N denote the submatrices of A and B of type r xr and 1 x r, which are obtained 
by simultaneous change of basis which puts A and B in block-diagonal form: 



A 



M 




B=\N ] . 



Simultaneous diagonalization is possible because of the imposed assumption BK = 0. 
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F Matrix theorems 

Here we explain some results for matrices that we have used in the main text. These results can 
be found in [11]. However, one of the results, the statement (c) bellow, is a new result, to the 
best of our knowledge. 

Theorem 1. Let A be a symmetric real matrix of type n x n and let R be its rank. Let us 
split A into blocks as 

' S N ' 

N'^ M 



A 



where S* is a J x J matrix, A^ is a J x r matrix, M is a r x r matrix and n = J + r. We will 
also assume that M is invertible. 

Let us construct the Schur complement (see [11]) S, which is a J x J matrix 

S = S- NM-^N^. 

Denote the rank of S as p. Then 

(a) R = r + p (Guttman rank additivity); 

(b) det A = detS'detM (Schur determinant formula); 

(c) if < /9 < J, then 



det Ma (det S4) 



det M det M^. 



(70) 



Here Ma and Mg are invertible R x R and p x p matrices, respectively. They are obtained 
by using orthogonal transformations which put A and S into a block-diagonal form 



A 





Ma 



S 





Mg 



' I 


-NM~^ ' 




S 


N ' 




I ■ 




' S 





I 




. ^^ 


M 




_ -M-^N^ I 




M 



The B4 matrix will be explicitly constructed in the proof below. 

Proof. We start from the Aitken block diagonalization formula [11] and from now on we use 
/ to denote a unit matrix of any size appropriate for its position in an equation 



(71) 



This equation can be verified by a direct multiplication of the left-hand side. Denoting the 
first matrix on the left as A, we can rewrite this identity in a compact form AAA^ = S (B M. 
The rank of the right-hand side is the sum of ranks of S and M, which amounts to p + r. Since 
the rank of A is equal to its dimension n, the total rank of the product on the left-hand side is 
equal to the rank of A, so we easily obtain 

R = r + p, 

which completes the proof of part (a). 

Next, we take the determinant of (71). Since A is block-triangular, its determinant is a 
product of determinants of blocks on the diagonal, so that we obtain det^l = 1. The left-hand 
side is thus the product of determinants, det^det Adet A''", and it is equal to det A because 
det^ = det74 = 1. On the right-hand side we have a block-diagonal matrix, so that its 
determinant is equal to det S det M. Hence, 

det A = det5detM, 
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which completes the proof of part (b). 

In order to prove (c), let O be a J x J orthogonal matrix which transforms S into a block- 
reduced form, 



OSO' 



0®M^. 



Since p / 0, matrix S has exactly p nonzero eigenvalues, which constitute M^, and since S is 
also real and symmetric, there will always exist an orthogonal matrix O that diagonalizes it. 
Given that the eigenvalues of M^ are nonzero, it is invertible. The zero-block is of type u x v, 
where v = J — p is the dimension of the null-space of S. By using O one can construct an 
orthogonal n x n matrix P = O ® I such that 



p ( 5 e M ) p^ = e M^ e M. 



(72) 



By using an analogous argument one can always construct an orthogonal n x n matrix Q 
such that 

(73) 



Q'AQ = OeMA, ^ A = Q{OeMA)Q. 



The zero block comes from the null-space of A. It is of the size n — R, which is also equal to i/, 
since n = J + r and R = r + p according to the part (a) of Theorem 1. 

Consider (71), and multiply it by P from the left and by P from the right, and use (72) 
and (73) to rewrite it in the form 



PAQ (0 e Ma) Q'A'P' = M^ M . 
Let us introduce the matrix B = PAQ and write it in the block form as 



(74) 



B 



Bi B2 
B3 B4 



where the blocks Bi, B2, -B3 and B4 are uxu^uxR^Rxv and Rx R matrices, respectively. 
Substituting this into the left-hand side of (74) yields 





Ma 



B 



T 



PAQ (0 Ma) Q' A'P' =B 

By comparing (75) to the right-hand side of (74), we obtain 



M^ 
M 



B2MABJ B2MABJ 
B4MABJ B^MaBJ 



(75) 



B2MABI B2MABI 
B^MaB^ B^MaBJ 



(76) 



Note that the zero-block of (76) is a i^ x i^ matrix, which is also the B2MAB2 block. We 
then read off the following equations 



B^MaBI = Ms.®M , 



B2MABI 
B2MABI 



0, 
0. 



(77) 

(78) 
(79) 



By taking the determinant of (77), we finally obtain 

det MA(det ^4)^ = det M det M^ . 
This establishes (70) and completes the proof of part (c) of the theorem. 



30 



Given that M, Mg and Ma are all invertible, we have deti?4 7^ which means that B/^ is 
also invertible. By multiplying (78) by (-BJ)~^M^ from the right, we obtain 

^2=0. 

The equation (79) now vanishes and does not provide any additional constraint. Therefore, the 
matrix B has the following form 



B = PAQ 



Bi 
B3 




Ba 



(80) 



End of proof. 

Remark 1. The A matrix from the main text has the form 



A 



S N 

N'^ M 





which differs from the one in Theorem 1 by an additional zero-block of size Q — n. However, 
these additional zeroes are integrated out before the Theorem 1 is applied, and they do not 
affect the statements of Theorem 1. 

Remark 2. The result (c) is a generalization of the result (b) to the case when A is a 
singular matrix. While the part (b) is in fact valid for singular matrices, it merely states that 
= and provides no information about nonsingular principal minors of A. The result (c) is 
more fine-grained, and provides precisely this nontrivial information about A. 

It was assumed in the part (c) that < p < J. li p = J then A is a regular matrix, and 
hence the result (b) can be used. If p = 0, then S" = 0, z^ = J, and instead of (77) we obtain 



B^MaBI = M , 



and consequently 



In this case we can set P 



for the matrix B. 



detMA(det54)2 
/ and obtain 



detM. 



B = AQ 



Bi 
B3 B4 



Remark 3. In the main text we use the results (b) and (c) to determine the leading A-order 
of the Schur complement S, knowing that A is of 0(1). However, it is necessary to show that 
-B4 is of 0(1) as well. In order to do this, note that 

detB = detPdet^detQ = ±1, 

since P and Q are unitary matrices. On the other hand, from (80) we know that det B = 
det Bi det B4, so that we have 

detBidetB4 = ±l. (81) 

Let us now assume that the blocks Bi and B4 are of order k and m in 1/A, respectively 
1 



B, 



C_ 

Jk 



O 



Afc+i 



Ba = 1^ + (—^ 

A™- V A"^+-'- 



k,m^O, C,Dr^O{l). 
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The numbers k and m cannot be negative since the whole B matrix must be of 0(1). Namely, 
the matrices P and Q are orthogonal, and consequently all their elements are bounded above 
by 1. Thus P and Q are of 0(1). The matrix A is also of 0(1), since A and consequently 
M,N,M~^ are all of the same order. Therefore, B = PAQ ~ 0(1). 
Since Bi \s a, v x u matrix and B4 is a R x R matrix, then 

deti?i = ^detO + o(^) , det B, = ^det D + O (^-^^ . (82) 

By substituting (82) back into (81) we obtain the consistency equation 

ku + mR = . 

Since both i/, i? > while k,m ^ 0, the only solution of this equation is k = m = 0. Therefore 

det54~54~0(l). 

In the case when i^ = the A matrix is regular and instead of the part (c) we use the part 
(b) of Theorem 1. However, the part (b) does not involve deti34, so that we need the above 
result only for u > 0. 
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